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ABSTRACT 

Measuring the albedo of an extrasolar planet provides insights into its atmo- 
spheric composition and its global thermal properties, including heat dissipation 
and weather patterns. Such a measurement requires very precise photometry of 
a transiting system sampling fully many phases of the secondary eclipse. Space- 
based optical photometry of the transiting system HD 209458 from the MOST 
(Microvariablity and Oscillations of STars) satellite, spanning 14 and 44 days in 
2004 and 2005 respectively, allows us to set a sensitive limit on the optical eclipse 
of the hot exosolar giant planet in this system. Our best fit to the observations 
yields a flux ratio of the planet and star of 7 ± 9 ppm (parts per million), which 
corresponds to a geometric albedo through the MOST bandpass (400-700 nm) 
of Ag = 0.038 ± 0.045. This gives a la upper limit of 0.08 for the geometric 
albedo and a 3a upper limit of 0.17. HD 209458b is significantly less reflective 
than Jupiter (for which Ag would be about 0.5). This low geometric albedo rules 
out the presence of bright reflective clouds in this exoplanet's atmosphere. We 
determine refined parameters for the star and exoplanet in the HD 209458 system 
based on a model fit to the MOST light curve. 

Subject headings: extrasolar planets: HD 209458; ultraprecise photometry 

1. Introduction 

Among the roughly 250 planets discovered to date around nearby star^, about two 
dozen are observed to transit. Multicolour optical photometry of the transits provides accu- 

"'^MOST is a Canadian Space Agency mission, operated jointly by Dynacon, Inc., and the Universities of 
Toronto and British Columbia, with assistance from the University of Vienna. 

^The Extrasolar Planets Encyclopedia: |http:/ / www.exoplanet.eu 
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rate measurements of the orbital incl ination, and hence th e mass and radius of the planet 
with respect to the parent s tar (e.g.. iKnutson et al.ll2007a[). Infrared meas urements of the 



eclipses in thermal emission (jPeming et al.l 2005a : Charbonneau et al. 2005 ) give brightness 



temperatures at specific wavelengths to constrain the atmospheric properties. Measurement 
of the optical eclipse of such a system provides a geometric albedo, which when combined 
with the thermal emission data, sensitively probes the exoplanetary atmosphere to test 
for high-altitude cl ouds (such as silicate condensates) and possibly even weather patterns 
feowe et al.]l2006ah . 



Models of the interior structure of giant exoplanets require appropriate boundary con- 
ditions - planetary mass, radius and surface temperature - which can be set by optical 
photometry of transiting hot extrasolar giant planets (EGPs). HD 209458b, which orbits its 
host star at a distance of only 0.047 AU every 3.5 days, the thermal equilibrium temperature 
for the dayside photosphere can be estimated as 



T,, = %{Rj2aY'^[f{l-AB)] 



11/4 



'1^ 



where the star and planet are considered blackbodies, T^, and R^, are the effective temperature 
and radius of the host star, the planet at distance a with a Bond albedo of As and / is a 
proxy for atmospheric thermal circulation. The Bond albedo, Ab-, is the fraction of total 
power incident on a body scattered back into space. 

The geometric albedo Ag is the ratio of brightness at zero phase angle compared to an 
idealized flat, fully reflecting and diffusively scattering disk. From Ag, one can derive the 
ratio of the stellar and planetary fluxes 



A., 



Rr 



(2) 



where Rp is the planetary radius and a, the separation between star and planet. 



If the composition of HD 209458b is similar to Jupiter, then the surface temperature 
of the exoplanet is the dominant variable driving the atmospheric chemistry, and hence, the 
albedo. This means that the albedo s hould be a function of both the star-planet separation, 
and the luminosity of the host star (IMarley et al.lll999l ). The relationship between albedo 
and surface temperature is a good tracer of chemical processes (such as the formation of 
molecules) in planetary atmospheres. 

HD 209458 is a n earby (47 pc), br i ght ( V = 7.65) Sun-like star (GOV) hosting a tran- 



siting hot EGP (e.g., ICody fc Sasselovl |2002| ) . Spitzer observations at 24 fim detected an 



eclipse at a depth o f 0.26%±0.046% whic h gives a brightness temperature (1130 ± 150 K) for 
the observed band (IDeming et al.ll2005al l. Estimating the equilibrium temperature is highly 
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model-dependent (ISeager et al.ll2005l: iFortnev et al 



2007bl ). Atmospheric models by lFortney et al.l (120051 ) predict an effective temperature of 1442 



20051 : iBarman et al.ll2005l : iBurrows et al. 



K, assuming An reradiation and solar metallicity. 

The brightness of the star and the short orbital period (~3.5 d) of the exoplanet made 
this planetary system an ideal candidate for photometry by the MOST (Microvariablity and 
Oscillations of STars) satellite. The objectives of the observations were: (1) to measure 
accurately time s of several consecutive t ransits to search for undetected Earth-mass planets 
in close orbits (IMiller-Ricci et al.l 120071 ): (2) to search for t ransits of nearly- Earth-radius 
planets at other periods in the nearly continuous data sets (ICroU et al.l 120071 ): (3) and to 
measure or set a limit on the optical eclipse of the exoplanet HD 209458b. 



The first results of our eclipse analysis (IRowe et al.l l2006al ). based on a MOST trial 
run of 14 days in August 2004 sampling 4 consecutive transits and 4 secondary eclipses, 
set upper limits on the eclipse depth of 53 ppm (parts per million) and on the geometric 
albedo of 0.25. We present here improved measurements of the eclipse based on the original 
photometry and a second 44-day MOST run in 2005, sampling an additional 12 consecutive 
transits and eclipses. The combined datasets give la upper limits on the eclipse depth of 16 
ppm and on the geometric albedo of 0.08. 



2. MOST Photometry 



The MOST satellite flWalker. Matthews et al.ll2003l : Matthews et al.ll2004f ) houses a 15- 
cm optical telescope feeding a CCD photometer through a single custom broadband (400 - 
700 nm) filter. From its 820-km-high polar Sun-synchronous orbit, MOST can monitor stars 
in its equatorial Continuous Viewing Zone for up to 8 weeks without interruption. 

Photometry of HD 209458 was obtained in two runs, during 14 - 30 August 2004 and 
1 August - 15 September 2005, for a total of 58 d ays of high-cadenc e photometry. The 
data were collected in MOST's Direct Imaging mode (IRowe et al.ll2006al ). where a defocused 
image of the star is recorded in a subraster of the Science CCD, with 1.5-s exposures sampled 
every 10 s. Two fainter stars were simultaneously observed, HD 209346 (A2, V = 8.33) and 
BD -|-18 4914 (AO, V = 10.6) in the same way. The latter was found to be a hybrid {6 
Scuti -|- 7 Dor) pulsator (IRowe et al.ll2006bl ). We rejected exposures with high cosmic ray 
fluxes which occur when MOST passes through the South Atlantic Anomaly (SAA), as 
well as data with background illumination values greater than 3000 ADU due to scattered 
Earthshine modulated with the satellite orbital period (see §2.2p . This resulted in data sets 
with 106,752 and 334,245 points and duty cycles of 81% and 89% respectively for the 2004 
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and 2005 runs. 



The light curve of the 2004 run is shown in Figure 4 of iRowe et al.l (l2006a( ). and the hght 
curve of the 2005 run is plotted in Figure [H This figure shows the raw photometry and the 
photometry after reduction as described by iRowe et al.l (l2006al ). The reduction procedures 
include the removal of the correlation between background light and instrumental stellar flux 
measurements, correction of low levels of crosstalk noise between the Science and Startracker 
CCDs (read out at different rates), and estimation of photometric errors as a function of 
the sky background level. The last two steps (described in detail below) were critical to 
i mprove the limit on th e optical albedo of the exoplanet an d also for the transit timing 
(IMiller-Ricci et al.l 120071 ) and transit search (ICroU et al.l 120071 ) conducted with this same set 
of photometry. 
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Fig. 1. — The 2005 MOST photometry of HD 209458. Top: The unprocessed light curve. 
Middle: The data after corrections for stray light correlations. Bottom: The final reduced 
light curve, after crosstalk corrections. The data have been averaged in 40-min bins for 
clarity. The large scatter in the first half of the light curve is because the stray light levels 
where comparitively much higher. 
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2.1. 



Crosstalk Corrections 



The effects of crosstalk noise on the photometry can be seen in the upper and middle 
panels of Figure [1] as small dips (~ 2-4 mmag) in the light curve at times of 16, 18, 20 and 
22.5 d. Crosstalk (or "video noise") occurs when the CCD contr oller electronic s interact 



MOST instrument, the Science and Startracker CCDs are located close beside each other 
in the camera focal plane, and are electronically isolated from one another. The aluminum 
structure of the spacecraft bus provides excellent grounding for the camera electronics, and 
crosstalk levels tested prior to launch were well below the original mission science require- 
ments. However, the transit and eclipse analyses of HD 209458 were not part of the original 
plan for MOST and these results are more sensitive to the intermittent crosstalk. 

Crosstalk happens when both the Science and Startracker CCDs are read out at the same 
time. The duration of frame transfer and readout for the Science CCD is about 1 s, while the 
Startracker has a readout time of ~0.1 s. A noise band approximately 10 pixels wide appears 
on the Science CCD (corresponding to the ratio of readout times) which can overlap the 
subraster containing the target star. This noise source does not obey Poisson statistics, which 
allows us to correct for it effectively. Without crosstalk, the noise for pixels used to estimate 
the sky background level in the subraster can be predicted from Poisson noise expected 
from the incident stray light plus read noise inherent to the detector electronics. This can 
be compared directly to the standard deviation of the pixel intensity values. Crosstalk 
introduces additional noise which changes the ratio of measured to predicted noise, which 
would be constant outside of crosstalk events. The correlation of background light and 
instrumental magnitude (see §2.21 below) does introduce changes in the noise ratio, but the 
time scale of this variation (the satellite orbital period of about 101.4 min) and the durations 
of crosstalk events (about 0.5 d) are so different that they are easily distinguished from each 
other. 

The positive signal added by crosstalk noise will traverse through the downloaded Sci- 
ence CCD subraster as the times of the independent CCD readouts move out of synch. The 
background level will be overestimated, but the PSF (Point Spread Function) fit to the stel- 
lar image is largely unaffected. Hence, the stellar flux is underestimated and the width of 
the crosstalk-affected pixel region moving across the CCD gives the characteristic sawtooth 
shape of the crosstalk artifacts seen in the upper two panels of Figure [TJ By correlating the 
background pixel scatter with the stellar photometry, the background level is rectified and 
the correct stellar flux value is recovered. The bottom of Figured] shows the 2005 photometry 
after this correction. In particular, notice the adjustment to the depth of the flfth transit 
during day 16, which coincided with a crosstalk event. 



with another electrical component due to incomplete 
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2.2. Photometric Error Estimation 

It is still unclear what causes the correlation of sky background to instrumental stellar 
magnitude readings in the data. To investigate this, we tested the validity of the predicted 
Poisson errors to the observed scatter, outside of times of crosstalk and times surrounding 
exoplanetary transits (about 5 hr every 3.5 d). We binned the data by sky background 
value, with an equal number of data points for each bin. We then calculated the expected 
scatter in each bin based on Poisson statistics and the measured standard deviation. Figure 
|2] shows the ratio of predicted to measured scatter as a function of the background level. 
At background levels of about 3000 ADU (Analog-to-Digital Units), the peak of the stellar 
PSF becomes saturated, since the 14-bit ADC (Analog-to-Digital Conversion) of the MOST 
CCD electronics has a saturation limit of 16 384 ADU. This causes the large discrepancy 
between theory and measurement seen at this threshold in the figure. At low background 
levels there is a linear increase up to about 1000 ADU, followed by a slower rising plateau 
until saturation. This change could be interpreted as a variation in CCD gain with signal. 

We have corrected the error estimates with a linear interpolation over the binned sky 
background values for each photometric measurement. This provides the correct weight for 
each measurement and the proper statistics for subsequent fitting to models of the light 
curve. 



2.3. The Phased Light Curve 

Figure [3] presents the 2004 and 2005 photometry plotted in phase with the orbital period 
of the exoplanet (see Table [1]). The top panel shows the unbinned data. The repetitive 
pattern seen in these data is due to the increased photometric scatter during intervals of 
highest stray light modulated with the MOST orbital period of about 101.4 min. (This 
appears correlated to the planets orbits because, by coincidence, the orbital frequencies of 
the MOST satellite and HD 209458b are in a near-harmonic ratio of 50:1.) The middle panel 
shows the data binned in 40-min intervals. The effects of stray light are less obvious here, 
since we use weighted averages based on the photometric errors described above in §2.21 The 
planetary transit is at phase 0.75; the eclipse (not visible) at phase 0.25. The bottom panel 
shows the data averaged in bins of width 0.04 cycle in phase (~3.4 hr) with la error bars of 
about 25 ppm (/zmag). 
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Fig. 2. — Corrected photometric error estimates as a function of background level. E^eas 
is the measured error based on the standard deviation of each bin and Eca\ is the expected 
error based on photon statistics. The sharp increase in the ratio of these two errors after 
3000 ADU is due to saturation of the detector. 



3. Modeling the Light Curve 

The light curve of a system with one transiting planet contains the following variations 
in phase with the exoplanet orbital period: (1) the transits themselves; (2) the eclipses; and 
(3) the changing flux from the planet as it goes through illumination phases during its orbit. 
There are also possible intrinsic variations in the star due to rotation al modulation , some of 
which may be synchronised with the close-in planet's orbit (see, e.g.. IWalker et al.l (120061 )). 



The Bond albedo can be written as, 



AB = Agq 



(3) 
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Fig. 3.— The reduced 2004 and 2005 MOST photometry of HD 209458 phased to the period 
of the exoplanet. Top: Unbinned data. The pattern outside of the phase of transit is due to 
the coincidental harmonic relationship (50:1) between the orbital frequencies of MOST and 
the exoplanet, so modulated scatter in the data due to scattered Earthshine is also in phase 
with the period of HD 209458b. Middle: The data averaged in 40-min bins. Bottom: The 
data in bins of width 0.04 cycle in phase. Note the different magnitude scales for the three 
plots. 



where Ag is the geometric albedo and q is the phase integral. From the definition of Ag 
in Equation [2] the geometric albedo must be observed at zero phase angle. The chemical 
composition of the atmosphere, including particle s ize, can strongly influence the phase 



integral, such as strong backscatter (for example, see lGreen et al.ll2003l ). 



First we place limits on the albedo based on comparison of observations adjacent to the 
eclipse. The planetary transit lasts for ~ 3.7 hours (and hence the eclipse). To measure the 
secondary eclipse for the planet disappearing behind the star, we average the photometry 
over phase bins of width 0.044 (~3.7 hours) centred 0.5 phase away from the transit (at 
phase 0.25 in Figure [3]) • We then find the average of the photometry from two adjacent bins 
with the same width. Our error is estimated by bootstrapping the means for each bin (see 
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§3.ip . This basic approach gives a flux ratio of the star and planet of 2.2 x 10 ^ ± 2.9 x 10 



If we make the assumption that the planet scatters as a Lambert sphere then q is 
simply a function of the phase of the planet from the observer's perspective. For a planet 
in a circular orbit like HD 209458b, we have approximated the reflected light variation as a 
sinusoid (adequate for the phase variations from a Lambertian sphere for the current set of 
observations). Our function to model the light curve, including transits, is 



/ 



MA {R 



* 5 "I ^•i^'i 



c„,t) + ^ [l + cos(^ + 0)] 



sin[i] 



1 + y-sin{i) 



(4) 



where MA is the normalized flux computed using the analytic expressions of iMandel fc Agol 



(I2OO2I ) for a transiting planet, F^, Fp, 



and Rp are the fluxes and radii of the star and 
planet, a and P are the semi-major axis and period of the planet's orbit, i is the inclination 
of the or bit relative to our line of sight, c„ are the nonlinear limb darkening parameters for 
the star (jClaretl l2000l ) and (f) is an arbitrary starting phase to coincide with the transit in 
the phase diagram. A schematic of this light curve model is shown in Figure |H The transit 
depth for HD 209458 extends to a relative flux of 0.98, far below the plotted scale. The 
amplitude of the flux variations due to phase changes of the planet is equal to the depth of 
the eclipse of the planet by the star. 

We adapt a Bayesian approach to our best-model-fit minimization so that we can in- 
corporate priors in our fits. The probability function that we wish to maximize is 



p{yi\i)-p{yn\i)p{D\yi, Vn, i) 

piD,I) 



(5) 



where yi to ?/„ are the 13 model parameters listed in Tabled! p(yi\I) to are the cor- 

responding priors, p{D\yi, yn, I) is the likelihood function and p{D, I) is a normalization 
factor. Having priors is important as single-band photometry gives few constraints on the 
orbital inclination of the exoplanet due to degeneracy of this parameter with the radii of 
the planet and star. Our ada pted priors for t he orb ital inclination, orbital period, mass and 
radius of the planet are from lKnutson et al.l (l2007al ). based on multi-band HST photometry 



with limb-darkening inf ormation to esti r nate t he inclination angle. We take our prior for the 
mass of the planet from iLaughlin et al.l (l2005l ). Limb darkening parameters for a non-linear 
model were determ ined from a Ku r ucz mo del representative of HD 209458 show in Figure O 
kindly provided by lKnutson et al.l (l2007a| ). scaled according to the MOST custom passband. 
For Sun-like host stars, the peak of the emitted spectral energy distribution is in the opti- 
cal range, so the MOST bandpass geometric albedo is a good approximation to the mean 
integrated value of the geometric albedo. 
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Priors are given by 



Exp 



2e? 



(6) 



where /ij is the prior value of i/i, with a la uncertainty of Cj, as hsted in Table [H We 
do not fit explicitly for the limb darkening parameters. In the case of fixed parameters, 
p{yi) = 6{yi — fii) and p{yi) = 1 when no prior information is available. The likelihood 
function is given by 



n ^ 

p{D\y,,...,yr„I) = Exp[-x'] J] -= 



1=1 



27rcr,- 



X 



E 

i=l 



{di - f{xi;yi,...,yn) f 



at 



where cxj is the photometric uncertainty as calculated in §2.2 



(7) 
(8) 



We find the maximum of equatio n [5] using a down hill simplex model based on the 
Amoeba routine in Numerical Recipes (jPress et al.lll992l ). We did not adopt any prior for 
the radius of the planet, as a direct comparison of MOST and HST photometry of the HD 
209458 transits, scaled to the san ae passband, pr o duces different depths of transit. Using 
the spectroscopic observations of iKnutson et al.l (l2007al ) and the MOST bandpass, HST 
photometry was rescaled to match MOST photometry. Figure [6] compares the averaged 
transit for the two sets of photometry. The MOST data indicate a deeper transit than the 
HST data. This may indicate a systematic error in the measured amplitude of thr transit 
from MOST data. It may also be related to the non-differential photometry for HST. Unlike 
MOST, HST cannot continuously observe HD 209458. Instead the complete transit much 
be pieced together from different orbits as the HST field-of-view will typically be occulted 
by the Earth during some portion of the orbit. This causes a thermal change in the optical 
telescope assembly which in turn changes the optical fo cus, producing the w ell documented 
breathing effect which is relfected in the photometry (IMakidon et al.l |2006| ) . Keeping the 
planetary radius as a free parameter, our best-fitted value for the model is given in column 
3 of Table [D 



3.1. Bootstrap Error Analysis 

To estimate uncertainties in our fitted parameters, we use a bootstrap technique. This 
involves randomly selecting data from the original time series and generating a new time 
series, with replacement. Replacement means that any data point can be chosen more 
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than once, but the total number of points is always the same, so some points will not be 
included in the new datasets. Each generated data set has a noise profile similar to the 
original time series and repeating the fitting procedure on a series of these randomised data 
sets produces a robust error distribution for our best-fitted parameter val ues. For other 



examp les and discussion of the bootstrap method, we refer the reader to ICameron et al. 
( I2OO6I ) and references therein. 

We performed ~20, 000 bootstrap iterations. Figure [7] shows the bootstrap results for 
some key parameters in our model fit. Table [T] gives la errors for all fitted parameters. 
These uncertainties were estimated by assuming a normal distribution and calculating the 
standard deviation of the bootstrap sample. 



4. An Upper Limit on the Albedo of HD 209458B 

Our best fit to the flux ratio of the planet and star {Fp/ F^) is 7 x 10~^ with a la upper 
limit of 1.6 X 10~^. Applying the best-fitted value to equation [2] gives the geometric albedo 
measured through the MOST filter: Amost = 3.8 ± 4.5%. 

The top panel of Figure [7| shows the bootstrap error analysis for Amost as a function 
of the planet radius Rp. The use of priors heavily constrains the system parameters for 
the planet (such as radius and mass) that are allowed. A glance at the errors stated in 
Table [1] reveals that these errors are unrealistically small. Therefore, we also tested our fits 
by removing the prior information. The derived errors ( also lis ted in Table [1]) are more 



in tune with other studies such as that of iKnutson et al.l (j2007al ). Photometry alone gives 



no information about the mass of the planet or star and our unconstrained fits that give 
uncertainities on the mass of the star are larger than the mass of planet from spectroscopy. 
The main subject in this paper is the flux ratio of reflected light from the planet compared 
to the host star and most of the other parameters are from fitting the shape of the transit, 
which applies to only a small fraction of the data. The flux ratio of the planet and star is 
largely independent of the transit fitting parameters. 

In the MOST bandpass (shown in Figure [5]) the Solar System giant planet geometric 
albedos are all greater than 0.4 and have a Bond albedo {Ab) greater than geometric albedo 
observed over all wavelengths (AgTox), If we assume that Ag ~ A^tot for HP 2094 58b and 



that 0.67 < AgTOTMs < 1 based on arguments presented in iRowe et al.l (l2006al ) we can 
estimate Ab < 0.12 at the la confidence level. 



Rayleigh scattering, Mie scattering and molecular absorption a re the dominant me cha- 
nisms that determine the reflected and emitted spectra of an EGP (IMarley et al.ll2006l ). In 
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the bluer portion of the optical spectra (wavelengths shorter than 600nm) Rayleigh scatter- 
ing in a clear atmosphere will reflect a large fraction of the stellar flux outwards. In the 
red portion (wavelengths greater than 600nm) photons will be absorbed deep in the atmo- 
sphere that will make the reflected spectrum relatively dark. The strong incident UV flux 
can produce a rich mixture of compounds from molecules producing a haze that can absorb 
incident U V photons and wil l darken the appearance of the planet in the blue part of the 
spectrum (IMarley et al.lll999l ). Jupiter and Saturn, which receive much less UV flux than 
HD 209458, contain unidentified species at the leve l of a few parts per 10^° that de crease the 



blue and green geometric albedos by a factor of 2 (jKarkoschkalll994l : iBurrows et al.ll2007al ). 
While eg uilibrium phot o chemistry in HD 209458b cannot produce long-chain hydrocarbon 
hazes (ILiang et al.ll2004j ) owing to the high temperatures, the importance of non-equilibrium 
photochemistry, or photochemistry of S and N, compounds has not yet been explored. 

In the absence of clouds, all hot Jupiter models do predict extremely low visible- 
wavelength geometric albedos, due to strong, broad absorption lines of neutral atomic Na 
and K. These atoms are known to suppres s the emitted visible- wavelength flux of brown 
dwarfs at similar atmospheric temperatures (lLiebertll200ll ). Some mode l atmospheres include 
self-co n sistent cloud format ion in a ID, complete cloud cover scenario (lAckerman fc Marley 
( l200ll ). ICooper et al.l (l2003l )). Nevertheless, the physics of cloud formation is a process that 
is not well understood or constrained both for particle size and densities which are governed 
by the competing effects of condensation and coagulation versus sedimentation (Marley et 
al. 1999). The low reflected fraction of incident radia tion readily rules out reflective clouds. 
New Spitzer observations at 3.6, 4.5, 5.8 and 8.0 /im (IKnutson et al.ll2007bl ) i ndicate a tem- 



pertur e inversion which requires an extra, unknown absorber at low pressures (IBurrows et al. 
2007bl ). Our low albedo limit means that, if the absorber is a cloud, then it must not be 
highly reflective. 

The efficiency of the planetary atmosphere to transport heat from the day to night side is 
commonly parametrized by the surface area over which the planet reradiates absorbed stellar 
flux. For a fully mixed isotropic atmosphere, the planet will reradiate over Air steradians, 
whereas, without circulation the planet only reradiates over 27r steradians. HD 20 9458b is 
expected to be tidally locked. Infrared measurements at 8 /im by I Cowan et al.l (120071 ) place 2 
a limits on the variations in IR flux with phase at 0.0015. Our low Bond albedo measurement 
requires that planet distribute at least 35% (at the 1 cr limit) of the absorbed stellar energy 
on the night side of the planet. Th is means that / from Equation [1] must be less than 2. This 
agrees with IBurrows et al.l (l2007bl ) who flnd models with circulation can give reasonable flts 
to Spitzer measurements. The atmospheric depth at which the bulk of the stellar radiation is 
absorbed determines whether advection or radiative transfer dominates energy t ransport and 



the efficiency with which a planet circulates heat to the nightside of the planet (jSeager et al. 
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20051 ). Figure [8] plots the Bond albedo versus the equilibrium temperature for different values 



of /. The horizontal line shows the 24 fim brightness temperature which gives a lower limit 
on Teg and the verticle line shows the 1 a upper limit on Ab- This constrains the equilibrium 
temperature to be 1400-?^ < Tgg < 16507^ depending on the efficiency of thermal ciruclation 
in the planet's atmosphere. 



5. Conclusions and the Future 

HD 209458b is less reflective than Solar System giant planets such as Jupiter. Our 
measurements place a la upper limit on Fp/F^ of 1.6 x 10~^. The inferred low Bond albedo 
< 0.12 rules out presence of highly reflective clouds in the atmosphere of HD 209458b and 
is also consistent with non-cloudy atmospheric models. It also constrains the equilibrium 
temperature to be between 1400 and 1650i^" The largest unknown in determining Teg is the 
efficiency with which the planet distributes heat from the dayside to nightside of the planet. 
Further infrared measurements at different phases of the planet orbit will place stronger 
constrains on this parameter. 

The MOST satellite re-observed the HD 209458 system in August - September 2007 for 
approximately 4 weeks. A new observing mode has been added that allows images to be 
stacked onboard the satellite before downloading. This dramatically increases the number 
of photon counts that can be recorded. Instead of obtaining a 1.5-second exposure every 
10 seconds one can obtain 1 exposure every 15 seconds composed of 15 stacked 1.5-second 
exposures. This has decreased scatter in the data by a factor of ~3. MOST will either detect 
the secondary eclipse or place a significance limit of a few percent on the albedo. Analysis 
of these data is underway. 
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Table 1. Fitted Parameters for HD 209458 Lightcurve. 
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Table 3. Table of parameters to describe the MOST HD 209458 light curve. 
Column (1) lists the 13 parameters used to describe the light curve. The priors 
in column (2) refer to constraints used to derive the best fit parameters listed in 
column (3). Column (4) gives the errors on the fitted parameters if no priors are 
used. Column (5) gives the units for each row in the table. 
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Fig. 4. — A model of the flux changes during the orbit of HD 209458b. The phase labeled 
A corresponds to the centre of the eclipse, when the flux is due to the star alone. Phase E 
corresponds to the centre of transit, where the normalized flux drops to about 0.98 (off the 
magnified scale of this plot). The sinusoidal curve (marked at phases B, C, D, F and H) 
is due to the changing illumination phase of the planet, as illustrated at the bottom of the 
figure. Phases Xi through indicate the ingresses and egresses of the eclipse and transit. 
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Fig. 5. — Total system throughput for the MOST sateUite optics and CCD detector is shown 
as the dashed hne. A Kurucz model with Teff = 6100 K and log g = 4.38 cgs, representitive 
of HD 209458a is also plotted. 
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Fig. 6. — Comparison of HST and MOST photometry of the transit of HD 209458b. The 
MOST data are open diamonds and the HST data (scaled to the MOST bandpass) as dots. 
The curve is our best-fitted model for the MOST light curve. 
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Fig. 7. — Bootstrap results for HD 209458 parameter fits to MOST photometry. Contours 
outline the 68% confidence region. Top: Planet radius vs. geometric albedo through the 
MOST passband. Middle: Constraints on the radii of the star and planet. Bottom: Con- 
straints on the masses of the star and planet. 
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Fig. 8. — The dayside Teq for HD209458b as a function of As for different values of / (see 
equation [1]) . The la hmit of Ab is shown as a vertical dashed line. 



